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' Abstract 



We investigate neutral pion electroproduction off protons in the framework of lieavy 
baryon chiral perturbation theory. The chiral expansion of the S-wave multipoles 
5^ i -Eo+ and I/o+ is carried out to three orders. There appear several undetermined low- 

energy constants. Three are taken from a recent study of the new TAPS threshold 
vr*^ photoproduction data, one is fixed from the proton Dirac form factor and the 
novel S-wave constants appearing at orders q^ and q^ are determined from a best 
fit (constrained by resonance exchange) to the recent NIKHEF and MAMI data at 
photon momentum transfer squared k"^ = —0.1 GeV^. The inclusion of a particular 
set of dimension five operators is forced upon the fit by a soft-pion theorem which 
severely constrains the momentum dependence of the longitudinal S-wave multipole 
Lo+ at order g^. We give predictions for lower photon virtualities for the various 
multipoles and differential cross sections. Further improvements are briefly touched 
upon. 
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1 Introduction 



Chiral perturbation theory (CHPT) allows to systematically investigate the 
consequences of the spontaneous and the explicit chiral symmetry breaking in 
QCD. With the advent of CW machines, pion production by real and virtual 
photons has become a major testing ground for predictions based on nucleon 
chiral dynamics. In particular, over the last years there has been considerable 
activity to precisely measure neutral pion electroproduction in the threshold 
region at various laboratories, in particular at NIKHEF (Amsterdam) and 
MAMI (Mainz). These measurements are performed at energies very close to 
threshold and typical photon momentum transfer squared of k'^ = —0.1 GeV^. 
After the pioneering measurement of Welch et al. [1] at NIKHEF, which con- 
centrated on the S-wave cross section, kinematically more complete differential 
cross sections for various photon polarizations e = 0.5 . . . 0.9 are now available 
based on the data of van den Brink et al. [2] (NIKHEF) and Distler et al. [3] 
(MAMI). Furthermore, data at lower \k'^\ have also been taken at MAMI [4]. 

In the review [5] we already considered charged and neutral pion electro- 
production off protons and neutrons in the framework of relativistic nucleon 
CHPT [6] to order (here, q denotes the small expansion parameter which 
can be an external momentum or meson mass and the order refers to the 
Lagrangian). Since then, the theoretical framework has been considerably re- 
fined based on the heavy mass expansion proposed in [7] and discussed in 
detail in the recent review [8] . In this framework and calculating the S-waves 
to order g^, we have already reexamined neutral pion photoproduction [9,10]. 
Novel P-wave low-energy theorems (LETs) have been found together with 
a good description of the new TAPS and SAL data for jp — » n^p [11] [12]. 
In particular, chiral loops are necessary to understand the small value of the 
electric dipole amplitude £'o+ at the n^p and 7r+n thresholds. Furthermore, 
the P-wave LETs have been shown to hold within 10% or better. While Pi 
can be inferred directly from the unpolarized differential cross section, for P2 
one has to measure polarization observables. A measurement of the photon 
polarization asymmetry is underway e.g. at MAMI. However, there are some- 
what model-dependent means to indirectly determine P2, see e.g. [13] and 
[14] . These LET predictions extended to the case of virtual photons [15] were 
also used in the analysis of the NIKHEF data [2]. The investigation presented 
here is thus a natural extension of the one for photoproduction presented in 
[9]. To orders q^ and q^, there appear two new counter terms with a priori 
undetermined low-energy constants (LECs). One can be fixed from a recent 
dispersion-theoretical fit to the nucleon electromagnetic form factors, i.e. from 
the radius of the proton Dirac form factor [16]. The other one will be deter- 
mined from a best fit to the differential cross sections at = —0.1 GeV^. 
As we will show, due to a soft-pion theorem, these fits are too constrained 
since they lead to a A;^-dependence of the longitudinal S-wave multipole Lo+ 
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only through Born and one-loop graphs. We therefore are forced to include 
the leading corrections to this soft-pion theorem in the counter terms which 
is formally of order in the Lagrangian. This in turn leads to a satisfactory 
fit of the existing data at k"^ — —0.1 GeV^ and we then make predictions for 
smaller values of where this dimension five operator is much less impor- 
tant. In addition, three /c^-indcpcndcnt low-energy constants are taken from 
the photoproduction calculation to order [10]. 

The paper is organized as follows. In section 2, we briefly summarize the perti- 
nent kinematics and deflne the various cross sections and structure functions. 
Wc also present the natural basis of S- and P-wave multipoles which appear 
in the transition matrix-element. In section 3, the chiral expansion of the S- 
wave multipoles Eq^ and Lo+ is given to three orders in small momenta and 
the flve combinations of P-wave multipoles to two orders. The corresponding 
low-energy theorems were already given in [15]. We also introduce an ap- 
proximation which facilitates the calculation of the energy-dependence of the 
S-wave multipoles. In section 4, we give the explicit expression for the various 
low-energy constants based on the resonance saturation principle. In section 
5, the best flts to the differential cross sections are presented together with a 
detailed analysis of the various multipoles at k'^ — —0.1 GeV^. In section 6, we 
give predictions for the new data taken at /c^ = —0.06 GeV^ and investigate 
in some detail the range of 0.04 . . . 0.06 GeV^. A short summary together 
with a discussion of possible improvements is given in section 7. The appendix 
contains some consideration about the S-wave LECs at order q"^. 



2 Formal aspects 

In this section, wc assemble all the necessary definitions for the cross sections, 
structure functions and multipoles. To be specific, consider the process 



where 7* is the virtual photon with k^ < 0. Denote by W the cm energy of the 
ttN system and its threshold value by Wo = + m, with = 134.97 MeV 
and m = 938.27 MeV the (neutral) pion and the proton mass, respectively. 
The following quantities are used 



7*W+P(Pi) 



7r°(?)+p(p2) , 



(1) 



e 



-1 
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(2) 



with q the pion cm momentum, ko the photon energy, e the photon polariza- 
tion, El the longitudinal photon polarization, ip the electron scattering angle, 
00 the pion energy in the tiN cm system and AW gives the invariant energy 
above threshold. As explained in some detail in [9], we will account for the 
pion mass difference Mt^+ — M^o in the loops but keeping one nucleon mass. 
The TT+n threshold is located at Uc — 140.11 MeV which will also be used as 
the charged pion mass. We introduce the ratio 

>0 , (3) 



and as a further dimensionless quantity of order one 



y varies between 0.93 and 1.12 for AW = ... 15 MeV. 

The unpolarized pion electroproduction triple differential cross section reads 

da aEf{W^ — m^) da da 

~ = Ty -jpr- , (5) 



dEfdQfdn^ in^Eimk^{e-l) dn^ dQ^ ' 



with Ty the conventional virtual photon flux factor, a = e^/iir = 1/137.036 
the fine structure constant and Ei/f is the laboratory energy of the incom- 
ing/outgoing electron. The differential cross section can be split into transverse 
(T), longitudinal (L), transverse-longitudinal {TL) and transverse-transverse 
(TT) terms. 



da 2Wq 



dVt^ W^ - m2 



Rt + eLRL + v2eL(l + e) cos0i?TL + ecos2(f) Rtt 



(6) 



where the Ri {I — T, L, TL, TT) are called the structure functions, cf) is the 
azimuthal angle between the scattering and the reaction plane, and the cor- 
responding polar angle, spanned by the photon and pion directions, is called 

cos 9 — q-k . (7) 
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One also uses the separated virtual photon cross sections 
dui 2Wq 



Rj, I ^T,L,TL,TT . (8) 



This completes the necessary definitions of kinematical quantities and cross 
sections. 

In what follows, we will consider the threshold region, i.e. small values of AW, 
typically AW < 15MeV and small photon virtualities, \k'^\ < 0.1 GeV^. In 
that case, the pion three momentum is small and it is therefore advantageous 
to perform a multipole decomposition. We confine ourselves to the S- and 
P-waves in what follows. Consequently, the current matrix element takes the 
form 



m 



AttW 



J = ia ^£'0+ + q - k Pi^ + ia ■ k q P2 + q x k P3 

+ ia • fcfc (^Lo+ -Eo+ + q-k (P4 - P5 - Pi - P2)) +ia-qkP5, 

(9) 

in terms of the two S-waves £'o+and Lo+ and five P-waves. We choose the 
following combinations of the more commonly used P-wave multipoles Pi+, 
Mi-t and Li±, where E,M,L stands for electric, magnetic and longitudinal, 
respectively and the ± refers to the total angular momentum of the pion- 
nucleon system, j = I ± 1/2 (with I the pion angular momentum). 



Pi = 3Ei+ + Mi+ - Mi_ , P2 = 3Ei+ - Mi+ + Mi_ , 

P3 = 2Mi+ + Mi_, P4 = 4Li+ + Li_, P5 = Li_-2Li+. (10) 

These are the combinations which appear naturally in the transition matrix 
element and lead to the most compact formulae for the various cross sections. 
All multipoles are, of course, functions of the pion energy and the photon 
four-momentum squared, like e.g. £^0+ = Eq^^uj, k"^). In what follows, we will 
drop these obvious arguments. The structure functions Ri can be expressed 
in terms of the multipoles as follows (in the S- and P-wave approximation) 



PT = |£o+ + cos^Pip + -sin2^(|P2p + iPap), 

Rl = \Lo+ + cos ^ P4P + sin^ 9 IP5P , 
Rj.^ = - sin ^ Rc[(Po+ + cos 9 Pi) P; + (Lo+ + cos 9 P4) P^] , 

RTT^lsinH{\P2\'-\P,\') . (11) 
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3 Chiral expansion of the multipoles 



To perform the calculations, we make use of the effective Goldstone boson- 
baryon Lagrangian. Our notation is identical to the one used in [9] and we 
discuss here only some additional terms. The effective Lagrangian takes the 
form (to one loop accuracy) 

rW , A-i) , ri^) , A^) , r{2) , Mi) /i 9^ 



where the chiral dimension (z) counts the number of derivatives and/or meson 
mass insertions. We will work out the S-wave multipoles to order and the 
P-waves to order as explained in some detail in ref . [9] . Here, the order refers 
to the various multipoles. The corresponding Lagrangians are one order higher 
since the photon polarization vector counts as order q. For the P-waves, we 
thus have to consider tree graphs with insertions from C^^^j^'^^ and one loop 
graphs with insertions solely from C^^li- For the S-wave multipoles, we have in 
addition to consider one loop diagrams with exactly one insertion from 
and tree graphs from C^^jf. In comparison to the photoproduction calculation, 
we have two additional terms with undetermined low-energy constants. One 
is from C^^]^ and is related to the Dirac form factor of the proton, FYik"^)-, 
and was already considered in the relativistic calculation in [5]. In the heavy 
fermion approach used here, the effect due to the finite charge radius of the 
proton manifests itself via 



to T 1 



+ Snp{\) , (13) 



with A the scale of dimensional regularization, the axial-vector coupling 
constant and = 93MeV the pion decay constant. The terms in the square 
brackets in Eq.(13) stem from the loops. The constant 6rip{X) can be fixed 
from the recent dispersion-theoretical fit to the nucleon electromagnetic form 
factors, < >^'^^p= 0.774 ± 0.008 fm^ [16]. The other counter term is from 
C^^lf and appears in the chiral expansion of the S-wave multipoles. We do not 
need the explicit form of this term in the effective Lagrangian in what follows 
and thus refrain from giving it here.Q 

Consider first the chiral expansion of the P-wave multipoles to order q^. This 
should be rather accurate for the large multipoles but is afflicted with some 
uncertainty for the small ones, compare e.g. the discussion by Bergstrom [13] 

^ As discussed later, we also need to take one particular term from C^jf since with 
the sole fc^-dependent counter term from >C^^ togther with the radius correction 
from C^^j^ one is not able to describe the existing data at rather large ~ 5M^. 
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for the photoproduction case. He finds a good description of the large multipole 
Mi_|_ — Mi_ but some sizeable deviation for the much smaller £'1+. As we will 
show later on, the existing electroproduction data are not yet accurate enough 
to pin down the small multipoles with great precision and we thus stick to the 
g'^ approximation. Q The chiral expansion of the P-wave multipoles thus takes 
the form 

p. _ ^Born ^ p.loop ^ ^ct ^ ^ = 1, . . . , 5 , (14) 



with pI°"^ = and Pf* = for i = 1,2,4,5. The Born terms include the 
contribution from the proton anomalous magnetic moment Kp, which to lowest 
order appears in the dimension two pion-nucleon Lagrangian. For the large 
multipoles Pi 2 3 we have the following Born (and counter) terms (in case of 
^3), 



pBorn ^ ^gjrNq(^^ . ) ^^^^3^ + ] ^ (15) 



pBorn_ ^9ttNQ 



^ (-(l + /t„)Vu;2_^2 + X (16) 

13cj^ - IScj^A;^ + 2u:'^Ml + Sfc^M^ + 5/€p(u;^ + k'^Ml - 2u:'^k'^)]^ , 



with (^TrAf = 13.4 the strong pion-nucleon coupling constant. We notice that 
the P-waves scale with the pion momentum q. The constant hp has been deter- 
mined from a best fit to the new TAPS data for 7p — >■ vr^p, hp = 13.0 GeV~^, 
a value which is well understood in terms of A and vector meson exchange. 
However, there are also new data from SAL for the same process. The analy- 
sis of these gives a somewhat higher total cross section above 7r~^n threshold 
[12]. To account for this, we also use the larger value of hp = 15.8 GeV~^ as 
determined in [9]. The small multipoles Pj^s have the following Born terms 



Born ^ ^ f 2 + 3^) y^T^ , (1^ 



P 



Born_ (^9nNq ^ 

SOTrm^ V uj"^ 



5 



'2x 

^jVa;2-P . (19) 



This should eventually be refined when more accurate data will become available. 
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All these Born terms are, of course, real. The one loop contribution to the Pi 
{i — 1, 2, 4, 5) takes the form 



p|oop _ 



\ly{y + pf 



^ fp ^p\y ^py ,y'^^ (p + %)(4p + p^ + 4y) ; 



(20) 



loop _ q 



pioop _ 



327r^m^ 



\ly{y + pf 



2 ,y 

IP +3 



py 



r 



+ \n--y 



3+ 8 3 + 12^^+ 3^ + 



ploop 
^4 



py 



2 1 



H{y,p) 



p ,y py y , - 

3 + 3"T-T + 3^/'-^ 



X [-P 



I r i> I 2\ 

y + ^ + y) 



{p + 2y)py 



H{y,p) 



ploop ^ eg^N^c q 1 

327r2m3 ^^(^ + p)3 

5 



3 + 3"T + r^^ 



l/(4p + p2 + 4^) 



with 



H{y, p) — arctan 



2x/y + p 



+ arcsm 



P + 2|/ 



H{y.p) 



V4p + p2 + 4|/ ' 



a; < a;^ 



(21) 



(22) 



(23) 



(24) 



Loop effects which renormalize physical quantities like e.g. the anomalous 
magnetic moment Kp or the pion-nucleon coupling are properly taken care 
of. For oj > ojc, these loop contributions become complex. This has to be 
accounted for by the substitutions 



y 



(25) 



and 



arcsin 



p + 2y 



TT 



V4pTp^T% 2 



+ i In 



p + 2|/ + 2^(y + p)(j/-l) 
V4p + p2 + 4y 



(26) 



For the imaginary parts of the P-wave multipoles, one can give a rather com- 
pact and accurate description based on the Fermi- Watson final state theo- 
rem. This serves as a consistency check and gives a handy and rather accurate 



8 



estimate about the size of the Im Pj. Of course, in the full calculations to 
be performed later, we are not using these approximate forms. Consider the 
rescattering diagram shown in Fig.l. The corresponding imaginary part can 
be cast into the product of a charged pion production process followed by a 
pion-nucleon charge exchange (CEX) reaction. Consequently 

Im {MfP = ql ■ Re(M)"^" , (27) 



where is a generic symbol for any one of the P-wave multipoles Ei_^., 
Mi±, Li± and J = 3/2 (1/2) for the 1+ (1—) multipoles, 02 j denotes the 
corresponding charge-exchange (CEX) scattering volume and q+ — 
To lowest order, these scattering volumina are given by 



„CEX _ „CEX _ 



Furthermore, the P-waves for 7p — > ir'^n can be easily evaluated in Born 
approximation from the pion pole diagram. 



V2F^ 



1 



2q- k 



{2 + p)Ml (2 + pyMt 



, (29) 



where the terms in the square brackets come from the expansion of the pion 
propagator. Evaluating the operator structure of Eq.(29) fixes the /^7^^". Com- 
bining these with CEX scattering volumina and the appropriate kincmatical 
factors, we have 



' ' ' 487r2F3M2 2 + p ' 

Im Pf' = Im Pf^ , Im pf^ = \ Im pf^ 



(30) 




71 



7U 



\ / 



p n p 

Fig. 1. Rescattering diagram. The solid circles subsume various subdiagrams. 
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This shows that P2 has the largest imaginary part. These approximations 
become very accurate as p approaches zero. Isospin breaking effects not taken 
into account in such a consideration are of the order (o;^ — M^o)^^^/Af^o ~ 2%. 

We now turn to the S-wave multipoles £"0+ and Lo+. Let <S be a generic symbol 
for either one of them. The chiral expansion carried out to order takes the 
form 

S = 5'^°''° + + + 5^* . (31) 



As already mentioned, S — S{u!,k^). However, from our photoproduction 
study [9], we know that in the threshold region, M^o < u < 150 MeV, 
the cj-dependence is mostly due to the unitarity cusp which is fixed by the 
Fermi-Watson final state theorem. Stated differently, this accounts for the 
dominant isospin breaking effect due to the sizeable pion mass difference, 
(M^+ — Mtto) /M^o = 18%. We will discuss the numerical implications of 
this approximation later on. Consequently 



and the chiral expansion, Eq.(31), is applied to S{k'^). This means that we can 
calculate S{k'^) at threshold in the isospin hmit. The Born terms are readily 
evaluated. 



^[{5 + 4:Kp)e -3{5 + 2Kp)M^]^ (33) 



+ 



8m2 



The order one-loop contributions are given by 



^^^2-jq3-loop ^ 



P 



+ 



1287r2mFj 



+ 

1 + p ' 2(l + p)3/2 

2 p 



arccos 



+ 



1 + p (l+p)3/2 



arccos 



-P 
2 + p 

-P 



2 + p 



(34) 



These are finite and can be checked against the result of the relativistic cal- 
culation presented in [5]. To get an idea about the uncertainties induced by 
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the approximation on the energy dependence discussed above, we have also 
evaluated to this order the full S{lj, k^). The difference of the full result and 
the approximation takes the form 



^0+ 



loop _ ^9ttN^ 




r Vy 


1287r2mF3 1 


ly + P 






P 

l + P 


loop _ ^9nN^^ 




r y'" 


1287r2mF3 1 


ly + P 


2 




P 



p-{2y + 3p)^l-y + 



4|/ + 4p + 
2Vy + P 



arccos ■ 



2 + p 



-4./1 



2-2^1^ + 
P 



P 



H{y,p) 



l+p (1 + 



arccos 



Vy + p 

4 / 



2+p 2+p 



(35) 



can 



We will come back to this later on. The order contributions to S{k'^) 
be grouped into the graphs which are proportional to qa and the ones propor- 
tional to g\. These are no longer finite and depend on the scale of dimensional 
regularization, A. The counter terms to be discussed below have an infinite 
piece to cancel these divergences and their scale-dependence is such that it 
cancels the one from the loops. The explicit form of these terms is (we first 
give the result of the g\ graphs) 



p/,2^ga3 _e^i^^//8 7p In 26 59p 5 



— 2(2-|-p) / dxJ 1 — + px{l — x) arccos , ^ - > , 
{ ^l + px(l-x)J 

^^^^ - 1287r3mF3l^3 6 ^ A 9+ 36 +''^6 2^ 



+ (5 + P) /i^l,V^+V^ 
^3 6^V P 2 



1 



+ j dx 





4(1 + p)x^ + (p^ - 2p - l)x - p(l + p)/2 - 2 
^Jl — x"^ + px{l — x) 



X 



X arccos . > , (37) 

^1 + px(l - x) ^ 

Similarly, one finds for the diagrams proportional to qa (these form an inde- 
pendent subset of gauge invariant Feynman diagrams) 
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} , 4(1 + p)x^ - (4 + 6p)x - 4 x 1 

+ / ax , — arcsm , > , 

{ ^1 - x2 + px{l -x) ^1 + px{l - x) J 

^ } x{2x-l) f ^ ^ , 7r(5 + 9pa;/2-4(l + p)x2^ 



i 1 — x"^ + px{l — x)\ ^ 



X^ + px{l -X)\ ^1 - X2 + p,x(l - X) 



2 + (p - 2)x - 2(1 + p)a;2 . x \] 

H 1 ^ V ^^^gj^ ^ J ^ . (39) 

Wl — + pa;(l — x) V 1 + P^(l — x)'^ ^ 

Notice that the last integrand contains an integrable singularity of the type 
— X for all p > 0. As a check, the last integral in Eq.(39) at p = gives 
TT^ — 37r + 9. The one-loop corrections to order within the approximation on 
the energy-dependence are thus determined. Finally, we turn to the counter 
terms at this order. First, there are the two terms from C^^lf which appeared 
already in the photoproduction case. We have refitted the sum oi + 02 since 
in Refs.[9,10] the full a;-dependence of Eq+ was considered. Within the same 
approximation used here, the value of ai + 02 is 7.85 GeV~^ as compared to 
^1 + ^2 = 6.60 GeV~^ from the full u;-dependence. We will also use this latter 
value as a measure for the theoretical uncertainty of our calculations. For the 
S-waves, we have a priori two new counter terms at order g^. 



E{kY = eM^{{a^ + a2)M^ - ask^} , 

L{kY = eM,{(ai + a2)M^ - a^^ + a^{Ml - e)} , (40) 

so that Lik"^) — E{k'^) ~ (1 + p). However, as proven in the appendix, in the 
soft pion limit one can show that 

aa + a4 = , (41) 



and thus there is only one LEG and furthermore a strong correlation between 
E{k'^Y^ and L{k'^Y^ at order q^. The low-energy constant 03 will be treated 
as a free parameter and pinned down by a fit to the available differential 
cross section data at k'^ = —0.1 GeV^. It turns out, however, that with a A;^- 
independent L'^*(A;^), i.e. with the /c^-dependence of L{k^) coming solely from 
the Born and loop graphs, one is not able to fit the existing data. We therefore 
have to include the first corrections to the soft-pion constraint Eq.(41) away 
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from the chiral limit. This induces terms of the type 



E^X,L^\^a,M'^e , (42) 



which are arising from terms in the Lagrangian and are thus of higher 
order. These are the minimal terms one has to take to be able to describe 
the data at fc^ = —0.1 GeV^. Of course, there are other counter terms at this 
order. These, however, merely amount to quark mass renormalizations of the 
already considered ^^-independent counter terms and will therefore be set to 
zero here. In the next section, we show how to estimate the pertinent low- 
energy constants from resonance exchange. A fully consistent 0{q^) calculation 
would also include two-loop graphs and one loop graphs with insertions from 
C)^lf (besides others). We take here the pragmatic approach and subsume all 
these effects in the effective couplings entering the higher order low-energy 
constants. In addition, the form factor correction due to the finite proton 
Dirac radius, cf. Eq.(13), contributes as follows, 

Eik')-'' = L{er' = - '^l^^f Sn,{X) , (43) 



and is completely fixed from the knowledge of < >Y [16], Srip{X = m) = 
6.60 ±0.29 GeV~^ using as determined from the Goldberger-Treiman rela- 
tion, qa = gTTNFn/m = 1.328.Q This contribution is always treated separately. 
We now turn to the estimate of the size of the resonance contributions to the 
{i = 1,2,3,4) and the P-wave low-energy constant bp. 



4 Resonance saturation of the low— energy constants 



In this section, we consider the resonance saturation hypothesis to estimate the 
numerical values of the various low-energy constants. This principle works very 
well in the meson sector, see e.g. refs. [17-19] and is also a good tool to estimate 
the LECs in the presence of baryons, see e.g. [8] and [20]. In the baryon 
sector, resonance saturation proceeds in two steps. First, the effective field 
theory contains mesonic (M) and baryonic (A^*) excitations chirally coupled 
to the Goldstone bosons and the nucleons. Considering these excitations as 
very heavy, but keeping the ratios of coupling constants to masses fixed, one 
produces a string of higher dimensional operators involving only pions and 
nucleons with coefficients given in terms of the resonance parameters. Second, 



This value is not very different from the one in the relativistic calculation, 
6rlf{X = IGeV) = 7.35 GeV-^ [5]. 
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one then performs the heavy mass expansion for the nucleons, i.e. 



£eff[t/, M, TV, TV*] ^ C,s[U, N] ^ C,s[U, N,] , (44) 



where denotes the velocity-dependent heavy nucleon field (in general, we 
suppress the subscript 'v'). 

First, consider t-channel vector meson exchange, in this case the coupling of 
the p° and the uj to the nucleon and the subsequent vector meson decay into 
the 77*^7* (a general discussion concerning the coupling of spin-1 fields to the 
pion-nucleon system has recently been given [21]). This leads to 



V j_ V _ V _ 4_ 9pN (1 + f^p) Gnp^ 9ujN (1 + f^u)) G-Kuj'y 

a,^a^-a^- -ay - ^^^^^ + 47rmM2 ' 

where the Vir^'-f couplings can be determined from the radiative widths TiV ^ 
7r°7), V = p^,u!. The tensor to vector coupling ratio ny is known to be large 
for the p, Kp ~ 6 . . .6.6, small for the a;, — —0.16 ± 0.01. The pN cou- 
pling constant is fairly well known whereas there is some sizeable uncertainty 
concerning Qi^^. To avoid these uncertainties, we use a simpler form based in 
part on the gauged Wess-Zumino action, using the KSFR relation and setting 
Kp = 6, = (for details, see [9]), 

V , V V V ^ uV ^ /Ac\ 

-^+-2=a,=-a,=^^^, bp-j^,, (46) 



leading to aY + ^ 2.67 GeV"^ and 6^ = 3.13 GeV"^ (see also [9]). 

Prom the s-channel (baryonic) excitations, the dominant one is A(1232) ex- 
change. We use the following ANn and AA^7 Lagrangians 



Can. = ^^^^A'^e^,(Z)a^7f • f^N + h.c. 
4m 



92 A 



Am? 

93 



A'^e^,(x)75F^^r3taA7v 



^^^A^e,,(X')^,(dxF''')T'^N + h.c. , (47) 
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with 



(48) 



and denote the relativistic spin-1/2 (Dirac) and spin-3/2 (Rarita- 
Schwinger) field, respectively. The third term in C/^n^ vanishes for real photons 
[22]. The off-shell parameters X, Y and Z are severely constrained from vrA^ 
scattering, the nucleon polarizabilities and the fit to the new TAPS data for 
neutral pion photoproduction, for details see Ref. [10]. Using Eqs.(47), one 
finds for the A contribution to the (to order in the Lagrangian) 



«i +«2 = IT. T^^i + mF(3 -2Z) 

247rm'^m^ I \ m/^ — m 

+m^{Y + Z + AYZ)j 
+^72(^(2X + 1)(1 - 2Z) + ^(1 +X + Z + 4XZ))} 

at = -at =at + at + ^^JyH^Z - l)[8^iF + g,{2X + 1)] 

+g3[m{2X' + V){2Z - 1) - 2mA(l + X' + Z + 4X'Z)]| . (49) 

As discussed in the previous section, we have to account for the first correction 
to the soft-pion theorem 03 + 04 = 0. This is entirely given by A-exchange 
and contributes to Lo+ as follows,^ 

L;;/ = -eM,/e2 (03 + 04) , (50) 



with 



(03 + 04)"^ 



ASnm^m'j^ 



gi 



2Y{4:Z - 1) - 



+ g2 



X 



(4Z - 1) + 



m/\ — m_ 
Sm'i — Am/^m + Zw? 



+g3 



Z -X' 



m 



4m(mA — m) 
{1 + X' + Z + AX'Z) 



(51) 



and similarly for i?( 



0+ 



^ We remark that this contribution to Lo+ is the one which really gives the first 
correction to the soft-pion theorem whereas the equivalent £"0+ term is a quark-mass 
renormalization of 03. 
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^ 487rm^m^ 



X 



9i 



-rri^iAYZ + AY + Z + -) - mm^i^YZ -AY + 2Z+-) 



4m/ 



+mi(28YZ -AY-7Z-7) ^UYZ + Y + Z + 1) 

m 



1 



(m — ttia)^ 



+92 



m{18YZ -7Y + 9Z - mA{30XZ - AX + 12Z + ^] 

+^(12XZ + 3X + 3Z + -) 

m 2 , 



m — TTT-A 



+5'3 



-%X'Z -AZ + AX' + 2 + ^^(1 + AX'Z + X' + Z)] ] (52) 

m J J 



A good check on the lengthy expressions Eqs. (51,52) (together with the anal- 
ogous terms of order not shown here) is that they fulfill the constraint 
L{k^)-E{k^)^{l + p). 



Furthermore, bp is the same as in the photoproduction case [9], 



^A _ 9-KNgi\/2 r 2m\ + rriAm - 
^ 12'Krin?m\ I 2(mA — rn) 

+ m(Y + Z + 2YZ)+mA(Y + Z + AYZ)'^ . (53) 

The couplings gi, g2 and the off-shell parameter X, Y, Z have been previously 
determined from neutral pion photoproduction, ttN scattering and the nucle- 
ons' electromagnetic polarizabilities. Nothing is known about the signs and 
magnitudes of and X'. These only show up in reactions with virtual pho- 
tons. In what follows, we will use these two parameters to obtain a best fit to 
the MAMI and NIKHEF data at k"^ = -0.1 GeV^. Since we do not constrain 
the magnitudes of (?3 and X'. this effectively amounts to a free fit with 03 ^ 04. 
We finally remark that g^ and X' also enter virtual Compton scattering and 
can eventually be pinned down within some broad ranges in the future. 



5 Fit to the existing data 



In this section, we show the combined fit to the NIKHEF (e = 0.67) [2] 
and the MAMI (e = 0.582 and 0.885) [3] data at k'^ = -0.1 GeV^ and the 
resulting multipoles. Note that the MAMI data are ^-integrated differential 
cross sections involving only Rt and Rl- In Figs.2a,b, the dashed lines give the 
best fit at order q'^, i.e. with the soft-pion constraint a^ + a^ = 0. Clearly, this 
does not describe the data. If one adds, however, the order counter terms 
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for the S-wave multipoles as discussed in the previous sections, one finds an 
acceptable fit (soUd fines) with a x^/dof of 2.296. In the resonance exchange 
picture, this amounts to 

53 = -125.7 ±6.2 , X' = -0.22 ± 0.09 . (54) 



Clearly, such a large value of the unknown third NA'y coupling constant 
indicates that one subsumes in its value other effects like e.g. from two- 
loop graphs. Alternatively, one can only take the q'^ counter terms and re- 
lax the constraint 03 = —04 in the fit. This leads to 03 = —1.37 GeV^^ and 
04 = —0.22 GeV^^, which are numbers of natural size. We remark that in a 
combined fit to the NIKHEF and the MAMI data, it is not possible to cor- 
rectly get the normalization of the two MAMI data sets for the two different 
values of the photon polarization e, compare Fig.2b. 



AW = 10 MeV 



AW = 12 MeV 



AW = 14 MeV 




AW = 4 MeV 



AW = 6 MeV 



AW = 8 MeV 




AW = 2 MeV 



0.5 



I I I I I I I I I I I I I I I I 




50 100 150 

B [degress] 



Fig. 2a: NIKHEF data. Open squares (circles): 4> = 180° (0°). 
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Fig. 2b: MAMI data. Left (right) panel: e = 0.885 (0.582). 

Tlie corresponding multipoles are sliown in Fig. 3 for AW = ... 15 MeV, all 
in units of 10^^/Af^+. For the S-waves, we give the real and the imaginary 
parts, indicated by the solid and dashed lines, respectively. We remark that 
Re Eq^ has changed sign as compared to the photoproduction case, it shows 
the typical cusp effect at the opening of the 7r+n threshold. In contrast. Re Lo+ 
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Fig. 3: Multipoles at k'^ — —0.1 GeV^. See text for notations. 
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is essentially energy-independent with a very small cusp. For the P-waves, we 
only show the real parts since the imaginary parts are very small. The large 
P-waves are all of the same size as shown in the lower left panel (Pi, ^2,-^3: 
solid, dashed, dash-dotted line, in order) whereas P4 and P5 are more than 
one magnitude smaller, see the lower right panel in Fig.3 (P4, P5: sohd and 
dashed line, respectively). 

In Fig. 4, we show the effect of the approximation on the a;-dependence. With- 
out refitting any parameters, we have added the correction Eq.(34) to the 
S-wave multipoles. The solid curves in Fig.4 refer to the approximation and 
the dash-dotted ones to the complete cj-dependence in the loops to order . 
Wc sec that the effect on £^0+ is small and somewhat more pronounced in 
Lo+. We remark that these curves should only be considered indicative since 
there is also a difference in the g^-loops, which we did not evaluate in detail. 
Furthermore, part of this effect would be absorbed in the values of the LECs 
which we used in the fit. 



Re EO* 



Re L0 + 




[ Q 

0.0 2.5 5.0 7.5 10.0 12.5 15.0 0.0 2.5 5.0 7.5 10.0 12.5 15.0 



AW [MeV] 



AW [MeV] 



Fig. 4: Approximation on the a;-dependence at order q^. 

Furthermore, we have also performed some fits with the larger value of hp — 
15.8 GcV~'^. It leads to a somewhat larger (in magnitude) value of = 
— 145.8±6.2 and X' = — 0.26±0.06, consistent with the value given in Eq.(54). 
The x^/dof is 2.945, i.e. worse than for the smaller value of bp = 13.0 GeV~^. 
To improve the fits with the larger bp, one thus would have to readjust the 
parameters in the photoproduction case. Clearly, this discrepancy deserves fur- 
ther experimental clarification. In what follows, we will always use the smaller 
value of bp. Wc would like to stress again, sec also [5], that for a good test of 
chiral dynamics, one needs data at lower photon virtualitics as witnessed by 
the large value we find for the coupling constant and the necessity of us- 
ing a particular dimension five operator. However, as decreases, the latter 
becomes less and less important. Finally, we note that the convergence in the 
S-wave loops and counter terms is slow, similar to the photoproduction case 
[9], whereas the Born terms at order can safely be neglected. 
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6 Predictions at lower photon virtualities 



Having fixed all parameters at k"^ = —0.1 GeV^, we can now make predictions 
for lower photon virtualities. We will concentrate here on the range of |/c^| 
between 0.04 and 0.06 GeV^ since at the upper end some data have been 
taken at MAMl and it is conceivable that in the near future one will not be 
able to go below 0.04 GeV^ since the pion has to leave the target.^ 

In Figs. 5 and 6, we show the real parts of the S-wave and the large P-wave 
muhipoles as a function of AW for P = -0.06, -0.05, -0.04 GeV^ (solid, 
dashed-dotted and dashed lines, in order). At threshold. Re Eq+ passes zero 
for ~ -0.04 GeV^. 
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Fig. 5: Predictions for the S-wave mulitpoles. For notations, see text 
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Fig. 6: Predictions for the P-wave mulitpoles. For notations, see text. 

For the new MAMl data, we show the various differential cross sections at 
A;2 = -0.06 GeV^ with e = 0.582 and AW = 2 and 8 MeV as indicated by the 



Predictions at lower photon virtualities can be supplied upon request. 
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solid and dashed-dotted lines, respectively. 
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Fig. 7: Predictions for the differential cross sections. For notations, see text. 
The S-wave cross section oq, defined as 

oo = l-Eo+l^ + eL l-^^o+l^ , (55) 



is shown in Fig. 8. The data of Welch et al. [1] (open squares) are taken at 
different values of e. This range is indicated by the dotted lines in Fig. 8. The 
solid line refers to e = 0.67, 
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Fig. 8: S-wave cross section. For notations, see text. 
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i.e. the kinematics of van den Brink et aL [2] (crosses). Within the sizeable 
uncertainty, the S-wave cross section as predicted by the calculation is 
consistent with the data and shows again the flattening with increasing 
as already found in the relativistic calculation [23] . 



7 Summciry and outlook 

In this paper, we have considered neutral pion electroproduction off protons in 
the framework of heavy baryon chiral perturbation theory. The main results 
can be summarized as follows: 

• Extending the successful photoproduction calculation [9] [10] to account 
for all A;^-dependent counter terms at order q'^ does not give a decent fit to 
the existing electroproduction data at the rather large k"^ — —0.1 GeV^ ~ 
— 5M^. This can be traced back to a particular soft-pion theorem which 
severely constrains the possible local dimension four contact terms. We 
are thus forced to take into account a particular dimension five operator 
which gives the first correction to the soft-pion theorem. With that term, 
the data can be fitted. 

• We have given predictions for the \k'^\ range between 0.04 and 0.06 GeV^, 
where new data have been and are going to be taken. At threshold, the 
real part of Eq^ changes sign at ~ —0.04 GeV^. The calculated S-wave 
cross section is in satisfactory agreement with existing determinations. 
More detailed predictions pertaining to the relevant kinematics of future 
experiments are available from the authors. 

Further improvements of the calculations presented here are: 

• The P-wave multipoles should be calculated to one order higher. This 
would give a more accurate description of the small P-waves and can 
further be used to tighten the predictions of the P-wave LETs [15]. 

• The full-energy dependence of the loops should be taken into account 
when more accurate data will become available. 

• A more systematic study of higher order effects than done here should 
be done to eventually pin down the A'"A7 parameters and X'. 

Clearly, when more and more accurate data at lower photon virtualities will be 
available, neutral pion electroprodcution offers a variety of tests of the chiral 
QCD dynamics. 
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A Counterterms at order q"^ 



Here, we wish to proof Eq.(41) which states that there are no polynomial 
S-wave counter terms with + up-to-and-including order g'^. For 

that, we follow appendix E of Ref.[5] and consider the soft-pion limit for 
the process 7*(A;) + p{pi) 7r°(g) +p(p2), i-e. for q^j, ^ (this includes the 
chiral limit 0). In this limit, only the nucleon pole graphs (including 

the electromagnetic form factors and, in particular, the radius counter term) 
remain, for all values of the photon virtuality fc^. The corresponding non-pole 
amplitude (denoted by an 'overbar') must thus vanish. 



+ {Bi + 2^4 - 2mBr) 



(A.l) 



where we use the notation of Ref.[5], i.e. P = {pi +p2)/2 (see section 2 of that 
review). In the soft pion limit, the following conditions arise 



B. 



Ai + B2 = , Si + 2^4 - 2mAe = . 



(A.2) 



In fact, the B^, i?i + 2B4 and Aq are individually zero since they are odd under 
crossing, i.e. proportional to (s — ti) ~ g • (pi +P2)- 

Modulo irrelevant prefactors and setting m — 1, the S-wave multipoles are 
related to the invariant amplitudes Ai and B^ as follows (we drop the subscript 
'0+' and refer to [5] for definitions) 



E^liAi+ iJ^As + I (a^^ - «^)^4 - i^Ae 
L -E^^ifi^ - iy)[-Ai -B2 + B1 + 2^4 - fiA^ -2Ae] , (A.3) 

with V k'^ p and \i ~ M^. To order q^^ the polynomial terms (in /c^) 
contributing to L, taking into account the crossing properties of the various 
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amplitudes, are solely given hy L — jiAi — v {Bi + 2Bi)/2, with 

1) Ai^ak^ and 2) Bi + 23^ = (3 {s - u)/2 = 2 (5 n . (A.4) 

Two cases are possible that could lead to 03 + 04 7^ 0, namely 

i) a or /3 , ii) a and /3 7^ with (3 . (A. 5) 

However, gauge invariance {J ■ k — Qi) (sec Eq.(2.1) of Ref.[5]) implies that 
B2 = —f3 k"^. Furthermore, the second condition from Eqs.(A.2) leads to B2 = 
—ak^ (because of 1) in Eq.(A.4)). Clearly, a and f3 must be equal and both 
possibilities in Eq.(A.5) are ruled out. This gives the desired result, 03+04 = 0. 

To summarize the argument, the soft pion theorem excludes a polynomial of 
the type Ai = const k"^ which would give rise to 03 + 04 7^ for the non-pole 
contributions. Only the case Ai = —B2 = const k'^ , B14 = const {s — u)/2 is 
allowed leading immediately to the constraint. Evidently, this constraint does 
not apply to the nucleon pole graphs, compare Eq.(43). 
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